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Abstract. We prove a monotonicity condition satisfied by the Erlang 
C formula when computed in the Halfin-Whitt regime. This property 
was recently conjectured in |Janssen et al.j [2011| . 



Recently, there has been a renewed interest in the Erlang C formula 

W ^ a ) = ^77^(E77 + ^-Vl > 

n\(l-p) k\ n\(l-p) J 

that gives the probability of waiting for an arriving customers to an M/M/n 
system whose traffic intensity is p = a/n. The value a is the offered load, 
that is the ratio between the arrival rate A and the service rate p. 

The Erlang C formula founds application in the dimensioning of large 
call servers, i.e. when the traffic intensity approaches the instability region 
(p — > 1) while the system keeps high its efficiency by accordingly increasing 
the number of servers (n — > oo). If the scaling is done in an appropriate 
way, the limiting system still shows a non degenerate behavior. 



This limit is known as the Halfin-Whitt regime, from Halfin and Whitt 



1981 , and it requires that the scaling is done via the square-root staffing 



principle, see also Borst et al. |2004| . The idea is to let the number of servers 



increase more than linearly with respect to the offered load a by adding a 
term proportional to the square root of a, i.e. 

(2) n(a) = a + /V« • 

In the scaling procedure the number of servers is allowed to take not 
integer values, hence it is useful to consider the extended version the Erlang 



C formula to the positive real numbers, see Jagers and van Doom 1 1986 

/ r°° nt 
(3) C(8,a)=\J (l + t y—e- at dt 
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that is valid for < a < s. In Halfin and Whitt 1981 , it is shown that 



(4) 



C*(/3) = lim C(a + /Va,a) 



1 m) 



-i 



with $ and (f) being respectively the distribution and the density functions 
of a standard Normal random variable. 

As observed in Janssen et al. [2011] , the value C*(/3) in Q can be used 
as a first approximation of C(a + fiy/a, a) for large a, and eventually as a 
dimensioning tool for large call centers as shown in Borst et al. [2004 , and 
supported by Janssen et al. J2009|. 



In particular in Janssen et al. 2009) , it has been conjectured that for any 
value of j5 the function C(a + j3y/a, a) is decreasing, that is the limit in (|4 



is 

approached from above and in particular the value of C*(/3) is a performance 
lower bound for any system with offered load a < 1. 

The purpose of this note is to prove that indeed this conjecture holds true. 
The main tool of the proof is in realizing that the function \/C{a + (3y/a, a) 
can be written in term of the moments of some special random variables 
that are proved to be stochastically ordered. It is interesting to note that 
even if the role of these random variables is crucial in the proof we were 
not able to give a physical interpretation to them. Probably if this could 
be done it would add additional insights to the structure of the function 
C(a + j3yfa,a) and in general to the asymptotic queueing system in the 
Halfin- Whitt regime. 

Theorem 1. The function C(a + (3y/a, a) is strictly decreasing in a for any 
fixed /3 > 0. In particular, for any a > 0, C{a + (3y/a,a) > C*(/3). 

Proof. Consider the function g(t, a) = at e~ at (1 + t) a_1 , with a, t > 0. It is 
non negative and J °° g(t, a) dt = 1 so we can look at it as a density function 
of a non negative random variable, say X a . 

Let Y a = (1 + X a )^, we can express the Erlang C function, see in 
terms of the inverse of the (3 moments of the random variables Y a , indeed 

C(a + (3Va-,a) = 1/E[(1 + X a f^] = 1/E[F/ ] , 

and we need to prove that the j3 moments, with j3 > fixed, are strictly 
increasing in a. Since the power function (-)P is increasing, it is enough to 
prove that the family of random variables {l^} a >o is increasingly stochasti- 
cally ordered. 

First we compute their density functions. Let y(x) = (1 + x)^, we have 

i . i i 

that the inverse x(y) = y^ — 1, with x > and y > 1 and x'(y) = ~^V^ , 

that is positive for any y > 1. It follows that the density function f(t,a) of 
Y a is given by 

f(y, a) = g(x(y), a) x'(y) = y^ a ~\y^ - l) e -«(^-D , 
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while its tail distribution has the following expression 

poo JL_ 

(5) F(y,a) = Pr{Y a >y}= / f(t, a) dt = e^^- 1 ) 

for y > 1. An easy check that Y^, with a > 0, is indeed a random variable 
can be done by differentiating equation ^ with respect to y and getting 
—f(y,a). Having that f(y,a) is non negative for y > 1 it follows that 
F(y,a) is a decreasing function in y, and the check is complete by noticing 
that F(y,a) — > as y — > oo and jP(1,o) = 1. To prove that the stochastic 
order holds it is sufficient to show that the function F(y, a) is increasing in 
a > 0. 

We can rewrite the tail distribution in the following way 



(6) F(y,a) = exp{(logy) 2 - (logy) 2 {f^fffi - 1)} 

log y log y 

= exp{(lo g2/ ) 2 M-^)} 
logy 

with h(x) = x + x 2 (l — e 1//x ). Since the exponential and the square root 
functions are increasing functions in their respective domains, we only need 
to show that h{x) is an increasing function for x > 0. Representing the 
exponential function by its Taylor series we have 

OO _ rt OO 9_n OO _ r) 

(7) h(x)=x + x ^-^ ) =x -^— = -Y.^ 

n=0 n=l n=0 v ' 

and the result follows because it is a sum of increasing functions. □ 



As a final remark, we mention that the Halfin-Whitt regime can be ob- 
tained in another way as well. That is, instead of writing the number of 
servers as a function of the offered load, like in ([2]), the latter is written as 
a function of the former in the following way 

(8) a(n) = n — (5y/n 

where in general (3 > 0, and the relation is valid for n > f3 2 . 

A natural question is to ask if the monotonicity property satisfied by 
C(s(a),a), a > is also valid for C(s,a(s)), with s > (3 2 . However the 
behaviour of the C(s,a(s)) looks more complicated as the following graphs 
show for two values of the parameter, f3 = 1/10 and (3 = 3. 
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